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Procedures  for  calculating  the  moisture  contents  of  1-  and  10-hour  timelag 
fuels  have  been  developed  based  on  theoretical  calculations  of  the  rate  of  mois- 
ture transport  in  wood.  The  1-hour  timelag  cal culation  is  superior  to  fine  fuel 
moisture  calculations  developed  previously  because  there  is  no  regional  bias, 
making  it  valid  over  a  wider  range  of  conditions,  and  because  it  separates  out 
the  effects  of  the  environmental  factors  of  temperature,  humidity,  and  solar 
radiation.  The  10-hour  timelag  cal  culation  produced  values  reasonably  consist- 
ent with  observations  obtained  from  1/2-inch  ponderosa  pine  fuel  sticks  exposed 
under  field  conditions. 


KEY  WORDS:   Forest  fuels,  forest  fire  hazard,  fuel  wood. 


o:.' 


Fuel  moisture  content  is  one  of  the  major 
variables  in  evaluating  fire  danger  and  pre- 
dicting fire  behavior.  Rothermel's^  develop- 
ment of  a  mathematical  model  for  predicting 
fire  spread  in  a  heterogeneous  fuel  and  its 
subsequent  adaptation  to  fire-danger  rating  re- 
quire moisture  inputs  for  more  than  one  class 
of  fuel.  Dead  fuels  have  been  classified  by 
their  moisture  timelag  for  fire-danger  rating .^"^ 
These  fuel  classes  are  the  1-,  10-,  and  100-hour 
classes,  which  correspond  roughly  to  cylindrical 

1/  Authors  are,  respectively.  Principal  Mete- 
orologist and  Research  Forester,  National  Fire- 
Danger  Rating  Project,  U.S.  Forest  Service,  Rocky 
Mountain  Forest  and  Range  Experiment  Station, 
Fort  Collins,  Colorado. 


fuels  less  than  1/4  inch  in -diameter ,"^/f  ;  to  1 
inch,  and  1  to  3  inches  in  diameter, 
respectively.  ^ 

In  the  past,  these  fuels  have  been  irbpSe- 
sented  by  physical  analogs;  basswood  slats  for 
fine  or  1-hour  timelag  fuels,  half -inch  ponderosa 
pine  fuel  sticks  for  the  intermediate  or  10-hour 
timelag  fuels,  and  2-inch  ponderosa  pine  dowels 
for  the  heavier  100-hour  timelag  fuels  (Gis- 

2/  Rothermel ,  R.  C.  A  mathematical  model  for 
fire  spread  predictions  in  wildland  fuels.  1971. 
(Unpubl i shed  report  on  file  at  N.  Forest  Fire  Lab., 
U.  S.    Dep.   Agr.,    Forest  Serv.,   Missoula,  Mont.) 

3/  Fosberg,  Michael  A.,  James  W.  Lancaster 
and  Mark  J.  Schroeder.  Dead  forest  fuels  charac- 
terization by  moisture  timelag.  (Manuscript  in 
preparation.) 
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borne  1936).  A  computational  system  is  needed, 
however,  because  field  use  of  physical  analogs 
is  not  always  feasible,  and  direct  measurements 
are  seldom  available  when  analyzing  records  for 
fire-danger  ratings.  Past  work  has  produced 
two  schemes  for  estimating  the  moisture  con- 
tents of  fuels  in  the  1-hour  group.  One  con- 
sists of  regression  equations  based  on  field 
measurements  of  basswood slats  (Storey  1965),^ 
the  other  consists  of  a  weighted  average  of  the 
equilibrium  moisture  value  for  wood  corre- 
sponding to  the  ambient  dry  bulb  temperature 
and  humidity,  and  the  observed  moisture  con- 
tent of  a  half -inch  ponderosa  pine  fuel  stick.  ^ 
In  general,  computational  schemes  have  not 
been  used  operationally  to  represent  the  10-hour 
timelag  fuel  moisture,  although  one  has  been 
developed  (Storey  1965).  The  100-hour  timelag 
fuel  has  generally  been  represented  by  a  buildup 
index  in  some  form.  Fosberg  (1971)  developed 
a  direct  computational  scheme  for  this  fuel. 

Previously  developed  computational  sys- 
tems were  generally  not  well  suited  for  universal 
application  because  the  derived  relationships 
were  based  on  limited  ranges  of  environmental 
conditions  characteristic  of  specific  geographic 
areas.  Thus,  for  nationwide  application,  the 
resultant  regression  equations  had  to  be  extrap- 
olated beyond  the  ranges  of  data.  This  limita- 
tion is  overcome  by  use  of  the  general  solution 
developed  by  Fosberg  et  al.  (1970)  and  Fosberg 
(1971). 

The  General  Theory 

The  basic  theoretical  solution  for  moisture 
gain  and  loss  is 


Am 


(1) 


where  the  relationship  between  change  in  mois- 
ture content,  5m=m.  ,  i-m.  ,  and  the  difference 

1  +  1  I 

between  the  surface  moisture  content  and  the 

initial  moisture  content,  Am=mb..   -m.  ,  is  a 

1  +  1  I 
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phase  tables.  National  F ire-Danger  Rating  System. 
54  p.  1966.  {Unpublished  report  on  file  at  Rocky 
Mt.  Forest  and  Range  Exp.  Sta.,  U.  S.  Dep.  Agr., 
Forest  Serv.,  Fort  Collins,  Colo.) 

5_/  U.  S.  Forest  Service.  Wildland  fire  danger 
rating,  n.d.,  n.p.  Pac.  Southwest  Forest  and  Range 
Exp.  Sta.,  Berkeley,  Calif. 


simple  exponential  function.  In  this  expression, 
is  the  moisture  content  at  the  end  of 
the  change  period,  i+1  ,  m^.  is  the  moisture  con- 
tent at  the  beginning,  and  mb^.^^  is  the  mois- 
ture content  of  the  surface  fibers  at  the  end. 
^  is  the  similarity  coefficient,  X  is  the  in- 
verse of  the  timelag,  and  t  is  the  time  period 
over  which  the  moisture  exchange  takes  place. 
The  similarity  coefficient  is  dependent  on  the 
product  Xt  .  C  is  derived  empirically  and  is 
used  to  insure  the  value  of6m/Am  which  is  con- 
sidered a  nondimensional  constant.  Stable 
solutions  exist  only  over  the  interval 
.05  <At<0  0.5.  This  implies  that,  for  1-hour 
timelag  fuels,  the  moisture  contents  may  be 
predicted  for  periods  of  one-half  hour  or  less. 
Since  we  are  interested  in  predictions  for  much 
longer  periods,  these  short-period  predictions 
must  be  assembled  sequentially.  To  do  this 
we  solve  equation  (1)  for  the  moisture  content 
at  the  end  of  a  time  step  1  +  1.  6t  is  the 
length  of  the  time  step;     t=  i  6 1  . 


5m  =  Am(l-^e"^^^) 


m. 


i+T  -  m.  =  (l-ce-^"^)(mb.^^  -  m.) 


m  =  (l-ce"^^^)(mb.^^  -  m . )  +  m. 


_  T,   X  -J- 

For  notational convenience,  let x=  1  -ce  ,thus 


m^.+T  =  x(mb.^^  -  m.)  +  m. 
m^.+l  =  xmb.^^  +  m.  (1-x) 


(2) 


To  solve  for  the  moisture  content  after  a 
number  of  time  steps,  i +0  ,  1,  2,--,  n-1,  n; 


and 


m^  =  xmb^.  +  m^(l-x) 


m^  =  xmb^  +  m^  (1-x) 

^-1  =>^-l  +m^.2(l-x) 

=  Xmb^  +  m^_i(l-x) 

where  the  initial  value  for  each  subsequent 
step  is  the  final  moisture  content  from  the 
previous  computation.  These  individual  solu- 
tions may  be  combined  to  give  a  solution  of 
the  general  form: 

i^=x[mb^+^?|]  (1-x)  ^*mb^_j  ]  +  ( 1  -x)  \  (3) 


m 
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Where  j  is  the  interval  of  summation;  j=n-i  . 

Equation  (3)  may  be  simplified  by  first 
considering  that  no  precipitation  occurs.  This 
assumption  may  at  first  seem  unduly  restrictive, 
but  the  effect  of  precipitation  can  be  added 
later. 

Since  drying  conditions  preceding  the  time 

for  which  the  fuel  moisture  is  to  be  evaluated 

are  cychc,  a  consideration  of  that  variation  must 

be  incorporated  in  the  computation.    This  is 

accomplished  by  defining  and  computing  a  clima- 

tological  coefficient  c    .which  is  the  ratio  of 

n-j 

the  equilibrium  moisture  content  at  the  end 
of  time  step  i,  me^  _ ^  and  the  equilibrium  moisture 

content  at  the  end  of  time  step  n,  me 


me 


n 


'n-J 


me. 


The  moisture  content  at  the  immediate  sur- 
face of  the  fuel  element  is  governed  by  the 
environment  and,  for  all  practical  purposes, 
is  the  equilibrium  moisture  content.  Except 
when  it  is  raining,  m^^  then  can  be  set  equal 
to  m^  .  Thus  ,  equation  (3)  becomes 


ran=x[me^+jl|  (l-x)-^mej^_j  ]  +  (! -x)  "m^ 


with 

me^  .  =  me^  c  . 
n-j        n  n-j 


m^=x[me^^me^  IJl]  (l -x) ^'c^.  .  ]  +    -x)  " 
m^=xme^  "X)  ^  c^.j  ]  +  ( 1  "x)  \  (4) 

and  the  term  (1  -x)*^  ni  becomes  very  small  and 
can  be  neglected.  The^precipitation  effects  are 
now  added  by  going  back  to  equation  (1)  and 
determining  a  moisture  change,  dmp  ,  due  to 
precipitation. 

dmp.+^  =  xmb.^^  +  m.(l-x) 

j6/  a  laboratory  experiment  showed  a  linear 
relationship  between  the  wetting  boundary  con- 
ditions mb  and  the  duration  of  wetting  td.  The 
constants  a  and  b  were  empirically  derived  from 
these  wetting  experiments  for  both  1/2-inch  and 
2-inch  ponderosa  pine  dowels. 

JJ  The  Kronecker  delta  is  a  mathematical 
notation  used  to  denote  a  situation  where  there 
are  two  or  more  possible  conditions  which  are 
mutually  exclusive,  i.e.,  valid  or  nonvalid,  ex- 
istent or  nonexistent.  In  this  case  it  denotes  rain 
or  no  rain. 


This  again  must  be  solved  in  series.  Thus 
wetting  becomes 
n-l 

dmp^  =  X.S  6.(l-x)'^  nib^_j  (5) 

6/ 

where  mb  is  the  wetting  boundary  condition, 

mb  =  a+b  t,.  where  t..     is  the  duration  of 
d  I  dl 

the  precipitation  during  time  step  1  .  The 
Kronecker  delta,-^     6^^    ,  is  ''0"  if  there  is 

no  precipitation  during  the  time  step  i  and  is 
if  there  is  precipitation.   The  dmp  given 
by  equation  (5)  is  added  to  the    m^  from  equa- 
tion (4)  to  give  the  total  moisture  content. 


Application  of  the  General  Theory 
to  Field  Problems 

The  solutions  provided  by  equations  (4) 
and  (5)  may  be  obtained  for  the  particular 
cases  of  the  1-  and  10-hour  timelag  fuels.  The 
solution  is  obtained  by  specifying  the  simi- 
larity coefficient,  C  ,    from  laboratory  studies, 

and  the  timelag,  X  \  the  time  increment,  6 1  , 
and  the  climatological  coefficient,     _ j  for  a 

1430  LST  observation.  For  the  1-hour  timelag 
fuel,  12  time  steps  of  V2-hour  duration  were 
used.  Thus,  ,  A=  1  , 6 1=0 . 5  ,  and  the  term 
X=0 .3935  .  The  climatological  coefficients  are 
averages  computed  from  data  taken  during  six 
general  observation  periods  of  the  O'Neill, 
Nebraska,  Great  Plains  Study  (Lettau  and  David- 
son 1957).  The  standard  weather  shelter  tem- 
perature and  humidity  readings  were  converted 
to  equilibrium  moisture  content  (U.  S.  Forest 
Service  1955).  These  coefficients  (fig.  1)  could 
then  be  substituted  into  equation  (4)  to  give 
the  resultant  equation 

m  =  1.0329  me  (6) 
for  a  midafternoon  observation. 

Since  the  equilibrium  moisture  content  de- 
pends only  on  temperature  and  humidity,  a 
table  readily  usable  in  the  field  may  be  con- 
structed (table  1).    Since  the  temperature  and 
humidity  are  measured  in  shelters  4V2  feet 
above  the  ground,  and  the  moisture  content  of 
the  fuel  depends  on  the  temperature  and  humid- 
ity of  the  air  immediately  in  contact  with  its 
surface,  the  shelter  readings  must  be  corrected 
to  account  for  the  temperature  and  moisture 
lapse  rates  between  the  levels  of  the  instru- 
ments and  the  surface  of  the  fuels  on  sunny 
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1.9  r 
1.8  - 
1.7  - 


1430      1330      1230       1130      1030      930  830 

LST 

Figure  1  .-The  dimatological  coefficients,  c^  j, 
for  a  1430  LST  observation.  The  data  points 
are  the  coefficients  calculated  from  the  bi- 
hourly  observations  for  six  observation  periods 
of  the  O'Neill,  Nebraska,  Great  Plains  study 
(Lettau  and  Davidson  1957). 


days.  The  temperature  and  moisture  lapse 
rates  depend  on  a  number  of  variables  not 
considered  because  of  the  complexity  that  would 
result.  They  are  windspeed,  aspect,  slope, 
stability,  and  the  radiation  absorption  and  emis- 
sion characteristics  of  the  underlying  surface. 
The  averages  of  values  found  in  the  literature 
result  in  corrections  of  15°  F.  increase  in  tem- 
perature and  a  3°  F.  increase  in  dew  point 
temperature  (Geiger  1957).  These  yield  an 
adjusted  relative  humidity  which  is  75  percent 
of  the  shelter  value. 

Increases  in  moisture  content  due  to  pre- 
cipitation are  difficult  to  compute  for  the  1-hour 
timelag  fuel  because,  first,  precipitation  duration 
must  be  known  to  the  nearest  half-hour,  and, 
second,  the  duration  must  be  applied  to  the 
particular  time  steps  in  which  it  occurred. 
The  second  difficulty  is  the  most  restrictive 
because  it  would  require  solution  of  the  series 
or  a  very  large  number  of  tables  to  provide 
for  all  contingencies.  Neither  of  these  are 
practical.  A  straightforward  alternative  is  to 
assume  that,  when  it  is  raining,  the  fuels 
are  at  fiber  saturation  or  30  percent  moisture 
content. 

Calculation  of  the  10-hour  timelag  fuel  mois- 
ture followed  the  same  procedures  except  that 
six  4-hour  time  steps  were  used.  Thus  6t=4  , 
A=0 . 1 ,  and  c,=  .98  ,  giving  a  value  of  x=  •  343 1 . 
The  dimatological  coefficients  were  determined 
from  the  same  source  and  in  the  same  manner 
as  for  the  1-hour  timelag  fuels.   This  gives  a 


 Table  1.— One-hour  timelag  fuel  moisture  (percent)  

State  of  weather^^  Relative  humidity  (percent)  

Code  0-1         Code  2-9       0     5     10      15     20     25      30     35     40     45      50     55     60     65     70     75     80     85     90  95 

  4.4.14.14.4.4.114,4.14.4.^^4.^4,100 

Temperature    Temperature    4     9     14     19     24     29     34     39     44     49     54     59     64     69     74     79     84     89     94  99 
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1/  In  recording  fire-weather  observational  data,  the  "state  of  weather"  code  0  indicates  a  clear  sky  and  1  indicates  5/lOs 
or  less  cloud  cover;  both  are  in  the  general  condition  "Sunny."    2  through  9  indicate  more  than  5/lOs  cloud  cover  and  various 
conditions  of  precipitation,  generally  included  here  under  "cloudy." 
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prediction  equation  for  the  moisture  content 
of  the  10-hour  timelag  fuels  of 


m  =  1 .2815  me 


(7) 


for  a  midafternoon  observation.  As  with  the 
1-hour  timelag  fuel,  shelter  readings  were  cor- 
rected to  account  for  temperature  and  humidity 
profiles  on  sunny  days. 

The  effect  of  precipitation  on  the  10-hour 
timelag  fuel  was  determined  by  considering 
precipitation  durations  of  1  to  4  hours  in  each 
of  the  six  periods  beginning  at  the  time  of  the 
observation  used  to  rate  the  day,  usually  in 
the  early  afternoon.  These  increases  in  mois- 
ture content  were  sufficiently  stable  to  allow 
the  24-hour  period  from  one  observation  to 
the  next  to  be  split  into  a  first  16-hour  period 
and  a  final  8-hour  period,  which  is  much  more 


practical  for  field  use  where  accurate  rainfall 
occurrence  records  are  not  easily  attainable. 

The  table  for  field  use  (table  2)  is  similar 
in  format  to  that  used  to  calculate  the  1-hour 
timelag  fuel  moisture.  If  precipitation  occurs, 
part  B  of  the  table  is  used,  and  the  correction 
for  precipitation  is  added  to  the  results  de- 
rived from  part  A. 

Certain  errors  result  when  the  computations 
are  made  with  data  from  the  tables  generated 
for  field  use.  The  final  term  in  equation  (4) 
is  not  always  negligible  for  the  10-hour  timelag 
fuel.  The  error  made  in  neglecting  this  term 
becomes  noticeable  only  when  the  moisture 
content  is  well  above  fiber  saturation,  a  situation 
that  is  not  important  in  fire-danger  rating,  how- 
ever, since  fuels  with  moisture  contents  above 
this  value  can  be  considered,  for  all  practical 
purposes,  to  be  fireproof.    Breaking  the  day 


Table  2. --Ten-hour  timelag  fuel  moisture  (percent! 
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Part  B-'^ 


Time 
precipitation 


Precipitation  duration  (hours' 


occurred         12       3      4       5       6  7 


9     10     11      12     13     14     15  16 


Basic  observation 

yesterday  to  0600       2       4       7       9     11      14     16      18     20     23     25     25+    25+    25+    25+  25+ 
today 

0600  today  to 

basic  observation       7     15     22     25+    25+    25+    25+    25+    -     -     -  -     -  - 

today 


ly  If  no  precipitation  has  occurred  since  basic  observation  time  yesterday,  today's  10-hour  timelag  fuel  moisture  is  read 
directly  from  this  part  of  the  table. 

2/  When  precipitation  occurs,  add  the  results  from  this  part  of  the  table  to  the  results  of  part  A.    The  sum  is  today's 
10-hour  timelag  fuel  moisture  value. 


5 


into  16-hour  and  8-hour  periods  instead  of  six 
4-hour  periods  for  computing  the  contributions 
of  rain  to  the  10-hour  timelag  fuel  moisture 
accounts  for  much  of  the  error  in  this  compu- 
tation. The  reason  is  that  the  magnitude  of 
the  precipitation  contribution  to  the  final  answer 
is  highly  dependent  on  when  the  precipitation 
occurs.  The  error  is  such  that  short-period 
rainfall  effects  are  overestimated  if  they  occur 
16  hours  or  longer  before  the  observation,  and 
are  underestimated  if  they  occur  within  4  hours 
of  observation  time. 


Evaluation  and  Comparison  of  the  Predictions 

The  1-hour  timelag  fuel  moisture  predictions 
have  not  been  compared  to  field  data.  Instead, 
they  are  compared  to  the  fine  fuel  moisture  cal- 
culations used  in  the  wildland  fire-danger  system 
used  in  California  and  to  the  fine  fuel  moisture 
calculations  developed  by  Storey  (1965)  and  used 
in  the  1964  version  of  the  National  Fire-Danger 
Rating  System.  These  comparisons  (fig.  2) 
show  that  predictions  developed  here  compare 
well  with  both  in  general,  but  that  under  condi- 
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Figure  2.-- 

Compari  son  between  predicted 
1-hour    timelag    fuel    mo  i  stures. 
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tions  of  low  relative  humidity,  they  are  more 
nearly  like  the  wildland  system  and  that  under 
periods  of  high  humidity,  they  are  more  like 
the  equations  developed  by  Storey.  This  is  a 
desirable  feature,  since  the  wildland  system 
was  developed  from  data  taken  in  arid  regions 
and  the  equations  developed  by  Storey  are  based 
on  data  collected  in  humid  regions. 

The  10-hour  timelag  fuel  moisture  calcu- 
lations were  compared  directly  with  the  half- 
inch  stick  observations  and  with  Storey's  pre- 
diction equation  as  modified  by  Schroeder  (1969) 
(fig.  3).  Comparison  with  the  Storey-Schroeder 
equation  shows  reasonable  agreement  except 


after  periods  of  precipitation.  This  is  to  be 
expected  since  the  Storey-Schroeder  equation 
does  not  consider  precipitation  duration.  Com- 
parison with  the  observed  half-inch  stick  mois- 
ture contents  is  also  reasonable  provided  one 
considers  (1)  the  limitations  introduced  by  the 
assumptions  made  in  predicting  the  contribution 
of  precipitation,  and  (2)  the  fact  that  the  half- 
inch  sticks  have  timelags  varying  from  12  to 
15  hours 

8,'  Personal  communication  with  William 
F i scher,  Northern  Forest  Fire  Laboratory,  Missoula, 
Montana. 


Figure  3.-- 

Compari  son  between  predicted 
and  observed  10-hour  timelag 
fuel  moistures. 
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Summary 
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from  1/2-inch  ponderosa  pine  fuel  moisture 
sticks. 


Fosberg,  Michael  A. 

1971.  Moisture  content  calculations  for 
the  100-hour  timelag  fuel  in  fire-danger 
rating.  USDA  Forest  Serv.,  Res.  Note 
RM-199,  7  p.  Rocky  Mt.  Forest  and 
Range  Exp.  Sta.,  Fort  Collins,  Colo. 

 ,  James  W.  Lancaster,  and  Mark  J. 

Schroeder. 

1970.   Fuel  moisture  response— drying  rela- 
tionships under  standard  and  field  con- 
ditions. Forest  Sci.  16:  121-128. 
Geiger,  Rudolf. 

1965.    The  climate  near  the  ground.  Re- 
vised ed.,  611  p.    Cambridge,  Mass.: 
Harvard  Univ.  Press. 
Gisborne,  H.  T. 

1936.  Measuring  fire  weather  and  forest 
inflammability.  U.  S.  Dep.  Agr.  Circ. 
398,  59  p. 

Lettau,  Heinz  H.,  and  Ben  Davidson,  [Ed.] 
1957.  Exploring  the  atmosphere's  first  mile. 
V.  2:  Site  description  and  data  tabu- 
lation,    p.  377-578.    N.Y.:  Pergamon 
Press. 
Schroeder,  Mark  J. 

1969.    Critical  fire  weather  patterns  in  the 
conterminous  United  States.  U.  S.  Dep. 
Commer.   ESSA  Tech.  Rep.  WB  8,  31  p. 
Storey,  T.  G. 

1965.  Estimating  the  fuel  moisture  content 
of  indicator  sticks  from  selected  weather 
variables.  U.S.  Forest  Serv.  Res.  Pap. 
PSW-26,  14  p.  Pac.  Southwest  Forest 
and  Range  Exp.  Sta.,  Berkeley,  Calif. 
U.  S.  Forest  Service.  Forest  Products  Laboratory. 
1955.  Wood  Handbook.  U.  S.  Dep.  Agr. 
Agr.  Handb.  72,  528  p. 


Agriculture— CSU,  Ft.  Collins 


8 


